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We report a magnetic force microscopy study of the magnetic domain evolution in the layered 
manganite La 2 - 2 ccSri+ 2 xMn 207 (with x = 0.32). This strongly correlated electron compound is 
known to exhibit a wide range of magnetic phases, including a recently uncovered biskyrmion phase. 

We observe a continuous transition from dendritic to stripe-like domains, followed by the formation of 
magnetic bubbles due to a field- and temperature dependent competition between in-plane and out- 
of-plane spin alignments. The magnetic bubble phase appears at comparable field- and temperature 
ranges as the biskyrmion phase, suggesting a close relation between both phases. Based on our 
real-space images we construct a temperature-field phase diagram for this composition. 

PACS numbers: 68.37.Rt, 75.60.Ch, 75.70.Kw 


I. INTRODUCTION 

Perovskite manganites with strongly coupled charge, 
spin, and lattice degrees of freedom provide a treasure 
chest of interesting and exotic physics.One of the 
most prominent phenomena studied in manganites is 
the colossal magnetoresistance (CMR).I^ In addition, a 
plethora of magnetic phases has been uncovered in the 
layered manganite La 2 - 2 ccSri+ 2 ccMn 2 07 , corresponding 
to the n = 2 member of the Ruddlesden-Popper se¬ 
ries (R, A)^+iMn^ 03 n+i {R and A being trivalent rare 
earth and divalent alkaline earth ions, respectively), such 
as antiferromagnetic metallic, ferromagnetic metallic, 
canted antiferromagnetic, insulating antiferromagnetic, 
and charge ordered.l^ The transitions among these states 
can be controlled by the doping level and are amplified 
in comparison to members of the n = 3 family due to the 
low dimensionality of the system, leading to a huge CMR 
effect. In particular, in the ferromagnetic phase, the easy 
axis can be tuned from out-of-plane (for 0.30 < x < 0.32) 
to in-plane (for 0.32 < x < 0.40).^^ Hence, the samples 
with X ^ 0.32 (“LSMO-032”) are close to an instabil¬ 
ity and a spin reorientation transition (SRT) can be in¬ 
duced easily by changing temperature and magnetic field. 
Several studies have been carried out to detail the mag¬ 
netic structure of LSMO close to the instability. These 
include low-field magnetic force microscopy (MEM) up 
to 240 zero-field Lorentz transmission electron mi¬ 
croscopy (TEM),I3 magneto-OTtical imaging,scanning 
tunneling microscopy (STM),and scanning Hall probe 
microscopy studies .^Intriguingly, a recent Lorentz TEM 
study revealed the formation of a biskyrmion lattice in 
thin LSMO samples under an applied magnetic field of 


about 0.3 T.EI However, until now there is no comprehen¬ 
sive report about the development of magnetic domains 
in a bulk single crystal of LSMO-032 exceeding magnetic 
fields of 240 GP 

Here, we report a detailed investigation of the domain 
evolution as a function of both temperature and mag¬ 
netic field in bulk LSMO-032 via MEM. We trace the 
SRT and find the formation of bubble domains in the 
vicinity of the SRT. In contrast to previous studies we 
directly observe the transition from dendritic to bubble 
domains as a function of the magnetic field. We also 
study the reversibility and unveil a weak hysteretic be¬ 
havior. Our temperature dependent MEM investigation 
evidences clearly the transition between bubbles and den¬ 
dritic domains. The MEM images reveal the gradual 
character of the SRT. We construct a magnetic field - 
temperature (H-T) phase diagram based on MEM im¬ 
ages obtained in applied magnetic fields of 0 T < /roH 
< 0.4 T and in a temperature range of 13.5 K < T < 70 
K. Our study serves as a counterpart to the recent inves¬ 
tigation of biskyrmions in a thin LSMO specimeiP^ and 
highlights the role of shape anisotropy in the formation 
of skyrmion phases in centrosymmetric compounds. 


II. EXPERIMENTAL DETAILS 

Single crystals of Lai.36Sri.64Mn207 were grown W 
the traveling-floating-zone method in an image furnace.^ 
Magnetization measurements were performed in a dc su¬ 
perconducting quantum interference device (Quantum 
Design). The dimensions of the sample under investiga¬ 
tion are 0.87(2) mm x 0.67(2) mm x 0.40(2) mm. The 
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single crystal was oriented via a four-circle diffractometer 
/ goniometer. 

MFM measurements were performed in a low temper¬ 
ature MFM system with a home-built MFM probe in¬ 
side a cryostat with a superconducting magnet in a field- 
and temperature range of 0 - 7 T and 4 - 300 K, re- 
spectively.l^ All experiments were carried out with com¬ 
mercially available tips (PPP-MFMR, Nanosensors). In 
MFM measurements, the magnetic force generated by 
magnetic domains is detected quantitatively in terms of 
the resonant frequency shift of the cantilever: a repulsive 
(attractive) magnetic force results in a positive (negative) 
frequency shift. 


III. EXPERIMENTAL RESULTS AND 
DISCUSSION 


We will first focus on the global magnetization re¬ 
sults before discussing our MFM images. Magnetization 
curves obtained at two different temperatures are shown 
in Figs. 1(a) (T = 90 K) and (b) (T = 20 K) with 
magnetic fields applied along the [100] direction (solid 
curves) and [001] direction (dashed curves). The satura¬ 
tion magnetization at T = 20 K and T = 90 K is 440 
emu/cm^ and 360 emu/cm^, respectively. No significant 
differences are found in curves between H // [100] and 
H // [110] (not shown), which can be expected consider¬ 
ing only a relatively small in-plane anisotropy.!^ In addi¬ 
tion, we do not find any evidence for hysteretic behavior 
in any field direction in the M-H curve. In order to 
directly compare the different magnetization directions, 
we use the corrected magnetic field Hcorr as the x-axis, 
which is obtained by subtracting the demagnetizing field 
Hd = D • M from the applied field H. The demagne¬ 
tizing factors D are based on the sample’s dimension. 
We find D = 0.22 for [100] and D = 0.49 for [001] .11^ 
From the magnetization curves we find an in-plane easy 
axis at 90 K and an easy axis along the crystallographic 
c axis at 20 K. Temperature dependent magnetization 
measurements in a weak applied field of /ioH = 5 mT 
yield the Curie temperature Tc = 105 K and a SRT 
from in-plane to out-of-plane at Tsrt = 63 K [see blue 
curve in Fig. 1(c)]. The absolute value of Tc together 
with the sharpness of the transition agrees well with pre¬ 
vious reports on LSMO-032 and underlines the sample 
quality in terms of doping homogeneity.!!^ The change in 
easy axis in samples with x ^ 0.32 originates from an 
additional contribution of the dipolar interaction (favor¬ 
ing an in-plane spin alignment) among Mn moments to 
the magnetocrystalline anisotropy constants Ki and K 2 
(fav oring an out-of-plane spin alignment for positive val¬ 
ues).!®^ To gain a deeper insight into the SRT process, 
we calculate the first and second order anisotropy con¬ 
stants {Ki and A 2 , respectively) from a fit applied to the 
hard-axis M-H curves at various temperatures (see Fig. 
1(c) and Ref.!!^ for details). For > 0 and A 2 > 0, 
the easy axis is along the c axis. However, for Ki < 0 
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FIG. 1. (Color online) Magnetization curves at (a) 90 K and 
(b) 20 K measured parallel (dashed lines) and perpendicular 
(solid lines) to the c axis, (c) Temperature dependence of 
the anisotropy constants Ki (empty black squares) and K 2 
(solid red triangles). Solid blue line: temperature dependent 
magnetization curve, with /xoH = 5 mT T c axis. 


and K 2 < —K\j2 the easy axis is within the ab plane. 
For < 0 and K 2 > —K\l2 the easy axis is along a 
cone centered around the c axis with an opening angle of 

d = sin“^ ^\/^7^ ^ case, above 80 K the easy 

axis is in-plane while below 40 K the easy axis is out- 
of-plane. The temperature evolution of the anisotropy is 
related to a change in orbital occupancy from at 

low temperatures to predominantly dx 2 _y 2 at high tem¬ 
peratures.!^ The small, positive Ki values below 40 K 
can be associated with a weak uniaxial anisotropy. Inter¬ 
estingly, a previously reported temperature dependence 
study of Ki in LSMO-032 found a strong increase with 
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lowering temperature, similar to our K 2 temperature be¬ 
havior P The fact that we ouly observe a weak, gradual 
iucrease of Ki as we decrease the temperature uuderliues 
the delicate balauce betweeu out-of-plaue aud iu-plaue 
spiu aligumeut iu our sample, as well as the seusitivity 
of LSMO-032 to tiuy compositioual differeuces. 

Iu order to elucidate the maguetic domaiu formatiou 
of the LSMO-032 siugle crystal iu au applied maguetic 
field, we image maguetic domaiu structures at T = 13.5 
K after zero-field cooliug (ZFC) through Tc- We start 
from the thermally demaguetized state aud apply a step¬ 
wise iucreasiug / decreasiug maguetic field iu the rauge 
of 0 - 0.36 T. A selectiou of these images is showu iu Figs. 
2(a) - (k). The greeu arrows deuote the measuremeut se- 
queuce. The domaius at zero-field exhibit clear deudritic 
structures with a uearly raudom orieutatiou (a weak ori- 
eutatiou prefereuce is observed, ruuuiug from the lower 
left to the upper right coruer), sigualiug ouly a weak 
iu-plaue auisotropy directed to 45° from the crystallo¬ 
graphic a aud b axes [see Fig. 2(a)]. This is iu accordauce 
with a previous maguetizatiou study, where a miuor iu- 
plaue auisotropy aloug the [110] directiou was revealed.^I^ 
These patterus result from au euergy miuimiziug pro¬ 
cess that compromises betweeu euergy reductiou from 
au overall decrease iu the surface maguetic-pole euergy 
aud the coucomitaut euergy iucrease from the creatiou of 
uew domaiu walls .1^ Similar domaiu patterus have beeu 
observed iu samples with a uuiaxial auisotropy.l^ Deu¬ 
dritic domaius remaiu at maguetic fields up to 0.18 T. 
For higher fields the maguetic laudscape displays au iu- 
homogeueous uucleatiou of bubble domaius (with a mag¬ 
uetizatiou opposite to the exterual field) iu coexisteuce 
with irregular stripe domaius [see bright spots iu Figs. 
2(b) aud (c)]. A crossover from the coexisteuce phase to 
isolated large bubbles (diameter 5 jam) occurs upou 
iucreasiug the maguetic field from 0.3 T to 0.355 T [see 
Figs. 2(c) - (e)]. At jj^oH > 0.355 T, a trausitiou to 
a uuiformly maguetized phase takes place [Fig. 2(f)]. 
This is supported by the maguetizatiou curve obtaiued 
at T = 13.5 K aud with H // c axis [Fig. 2(1)], where a 
sat ur at iou is reached arouud 0.35 T. Upou decreasiug the 
maguetic field agaiu, we observe the reappearauce of bub¬ 
bles arouud 0.35 T followed by au iuhomogeueous but re¬ 
versible trausitiou from bubble to stripe domaius arouud 
0.18 T [see Figs. 2(g) - (j)]. Note, that the directiou of 
the weak iu-plaue auisotropy is preserved, as evideuced, 
e.g., iu Fig. 2(h). Decreasiug the maguetic field further 
leads to the returu of deudritic structures. Here, we poiut 
out that although qualitatively the deudritic structures 
iu Fig. 2(a) aud (k) are the same, the local patterus 
are differeut. A similar behavior is observed at mag¬ 
uetic fields arouud 0.33 - 0.35 T, where bubble domaius 
vauish (upou iucreasiug fields) aud uucleate (upou de¬ 
creasiug fields) at differeut locatious. These observatious 
highlight the abseuce of chemical gradieuts but a homo- 
geueous stoichiometry throughout the sample, suggest- 
iug that the domaiu structure aud the uucleatiou sites 
are uot related to local defects. Ou the other haud, some 


features that cau be associated with metastable bubbles 
remaiu at zero fields [see arrows iu Fig. 2(k)]. Iu Fig. 
2(m) histograms of the 0 T MFM images before (solid 
red bars) aud after (empty blue bars) the field-cychug 
are compared. Both histograms show a similar distribu- 
tiou of uegative aud positive Auj values ceutered arouud 
Auj = 0, which agrees very well with the abseuce of auy 
hysteretic behavior iu the M-H curve. At the same time, 
the distributiou is cousiderably uarrower after the field- 
cooliug, which cau be directly observed iu the reduced 
coutrast of Fig. 2(k) with respect to Fig. 2(a). We cau 
therefore attest a reductiou iu the maguitude of the out- 
of-plaue maguetic momeut aud heuce a small maguetic 
field hysteresis iu LSMO-032. This was uot picked up 
iu previous maguetizatiou measuremeuts aud remaius au 
opeu issue.l^ The same sequeuce of images was also takeu 
at T = 50 K (uot showu), where bubble structures ap¬ 
peared iu the same maguetic field rauge. Regardless of 
the temperature aud the maguetic field, uo loug-rauge or¬ 
dered bubble-lattice structures have beeu observed. This 
is coutrasted by the Loreutz TEM study, where a well- 
arrauged lattice of biskyrmious has beeu reported iu thiu 
LSMO specimeus with x = 0.315.^ 

Iu the followiug, we iuvestigate the temperature depeu- 
deuce of the maguetic domaiu structures. The sample 
was field-cooled through Tc with a fixed maguetic field 
of jioH = 0.23 T, which left the sample with the coex¬ 
isteuce of stripes aud bubbles below Tsrt- The temper¬ 
ature depeudeuce of the domaiu structure of LSMO-032 
iu zero- aud weak maguetic fields {jHoH = 24 mT) has 
beeu iuvestigated previously via MFM, without observ- 
iug auy sigus of bubble domaius.!^ lustead, a chauge from 
iu-plaue (with clear domaiu walls) to deudritic structures 
arouud Tsrt = 85 K was observed. This is cousiderably 
higher thau our Tsrt of about 63 K at zero fields aud 
further highlights the delicate depeudeuce of the sam¬ 
ple’s maguetic properties ou miuute differeuces iu com- 
positiou arouud x = 0.32. The deudriti c do maius show 
drastic chauges betweeu 63 K aud 45 K.^EH This evolu- 
tiou cau be associated with a ripple state that evolves 
liuearly duriug the spiu reorieutatiou.^l The ripple state 
is related to the layered crystal structure that euables 
dipolar iuteractio us b etweeu Mu-ious (iu-plaue) aud the 
orbital trausitiou.^^^^ Figures 3(a) - (g) show the evolu- 
tiou of the maguetic domaius upou cooliug the sample 
from T = 66 to 34 K, aud Fig. 3(h) plots the tempera¬ 
ture depeudeuce of the maguetizatiou with H // [001] as 
well as T [001], aud fi^H = 0.23 T. At this field, Tsrt 
is reduced from ^ 63 K at = 5 mT [see Fig. 1(c), 
blue hue] to ^ 50 K. Ou the other haud, our MFM im¬ 
ages reveal that the spiu reorieutatiou process is rather 
gradual aud occurs over a wide temperature rauge. For 
temperatures betweeu 80 K aud 64 K, the sample hosts 
iu-plaue domaiuJ^ or couically aligued spius betweeu the 
ab plaue aud the c axis,l^ heuce produciug MFM images of 
weak coutrast. Weak domaiu structures emerge arouud 
62 K, aud from 56 K to 42 K the coutrast iucreases dra¬ 
matically, sigualiug au iucrease iu out-of-plaue momeuts. 
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FIG. 2. (Color online) (a) - (k) Magnetic field evolution of the domain structures at T = 13.5 K. Magnetic fields are indicated 
in the figures. The scale bar corresponds to 5 /rm. The crystallographic orientation for all images is indicated in (f). (1) 
Magnetic hysteresis loop at T = 13.5 K with H // c axis. The following phases are indicated: (i) dendritic domains, (ii) 
irregular stripes and bubbles, (iii) isolated bubbles, and (iv) saturation, (m) Histograms of the 0 T MFM images before (solid 
red) and after (empty blue) field cycling. 
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FIG. 3. (Color online) (a) - (g) 
Evolution of magnetic domain struc¬ 
tures in LSMO-032 with tempera¬ 
tures ranging from 66 K to 34 K 
and fioH = 0.23 T along the c axis. 
The scale bar corresponds to 5 /xm. 
The crystallographic orientation for 
all images is indicated in (a), (h) 

Temperature dependence of the mag¬ 
netization with iiqH — 0.23 T ap¬ 
plied perpendicular (blue circles) and 
parallel (red squares) to the c axis. 


Around 60 K, stripe domains and bubbles begin to co¬ 
exist [see Fig. 3(d) - (g)]. Here, two regimes are ob¬ 
served: First, between 60 K and 45 K [see Fig. 3(e)] 
disordered structures coexist with stripes and bubble do¬ 
mains. Subsequently, below 45 K [see Figs. 3(f) - (g)], 
the images show disordered bubble structures with fea¬ 
tures similar to those observed at 13.5 K at comparable 
magnetic fields. As evidenced in Fig. 2, a coexistence of 
stripes and bubbles even remains at intermediate fields 
(^0.2 - 0.3 T) down to low temperatures. Interestingly, 
the nucleation of bubbles segmenting stripe domains can 
be directly observed [see black dashed boxes in Figs. 3(d) 
and (e)]. 

In Fig. 4 we summarize our MFM and magnetiza¬ 
tion results in a schematic field - temperature (H-T) 
phase diagram of bulk LSMO-032. A similar evolu¬ 
tion from stripes to magnetic bubble domains can be 
expected in other bulk systems governed by a compe¬ 
tition between ferromagnetic exchange and long-range 
dipole-dipole interaction, such as M-type hexaferrite,^ 
as well as in artificially fabricated multi-layered sys¬ 
tems with opposing magnetic anisotropies, such as 


0.5 

0.4 

P 0.3 
X 

o 

0.2 

0.1 

0.0 



0 20 40 60 80 


Temperature (K) 


FIG. 4. (Golor online) Schematic phase diagram of a LSMO- 
032 single crystal with the magnetic field applied along the 
c axis. The white squares correspond to the SRT obtained 
from magnetization curves. The geometrical demagnetization 
factor along this direction is D ^ 0.49. 
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(Fe/Ni)/Cu/Ni/Cu(001)P^In the latter, the sandwiched 
out-of-plane magnetized Ni-layer competes with the in¬ 
plane magnetization of the top (Fe/Ni) layer via inter¬ 
layer coupling. A transition between stripes and bubble 
domains can be achieved by varying the thickness of the 
separating Cu-layer. In contrast, LSMO-032 is located 
very close to an instability, which allows for an easy do¬ 
main tuning via temperature and considerably low mag¬ 
netic fields without changing the composition. 

The composition of our specimen is also close to the 
one of the thin sample, in which a lattice of regularly ar¬ 
ranged biskyrmions with ^ 100 nm diameter and the 
topological charge Nsk = 2 was reported (with x = 
0.315).^ Although we cannot extract any information 
about the spin structure at the nanoscale from MFM, 
we classify our observed bubbles as topologically triv¬ 
ial magnetic bubbles with Nsk = 0- This assumption 
is supported by the fact that the bubbles are randomly 
arranged, appear to have varying diameters and shapes, 
and are significantly larger than the biskyrmions. Hence, 
the bubble formation is not (directly) related to the un¬ 
derlying crystallographic structure or to a delicate bal¬ 
ance of competing interactions and anisotropies. On the 
other hand, the similarity in composition between the 
two samples suggests the possibility of observing a (bi- 
)skyrmion phase in our sample in the presence of an 
additional shape anisotropy. Indeed, Yu, et ah, report 
a decisive dependence of the biskyrmion diameter on 
the sample thickness (i.e., on the shape anisotropy) 
This possibility opens the door of observing directly a 
crossover from topologically trivial magnetic bubble do¬ 
mains (in the bulk crystal) to topologically non-trivial 


skyrmions (in the thin sample limit) in the bilayer man- 
ganite LSMO-032 with a thickness gradient. 


IV. SUMMARY 

We presented a detailed investigation of the magnetic 
domain evolution in a LSMO-032 bulk sample in real 
space using MFM. Our field-dependent measurements 
uncovered the formation of bubbles from dendritic do¬ 
mains, with reversible domain structures but hysteretic 
magnetic magnitudes. In addition, our temperature de¬ 
pendent study shows the possibility of generating bubble 
domains through the field-cooling, as well as the grad¬ 
ual character of the SRT. Our results pave the way for 
future studies to reveal the relation between magnetic 
bubbles and skyrmions in materials with a high uniaxial 
anisotropy. 
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